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PRELIMINARY NEUTRONIC ANALYSIS OF A CAVITY TEST REACTOR 
by Charles L. Whitmarsh, Jr. 

Lewis Research Center 

SUMMARY 

A neutronic analysis was performed on a cavity test reactor to be used to test the 
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gas-core nuclear rocket concept. This reactor would provide a thermal flux of 4. 1x10 
neutrons per square centimeter per second to a 60.96-centimeter-diameter centrally 
located spherical test cavity in order to produce 2. 7 megawatts in the 375 grams of en- 
riched uranium plasma fuel contained there. A major constraint imposed on this study 
was that fuel elements of the type used in the Materials Testing Reactor (MTR) be used 
for the driver fuel. 

The analysis was directed toward minimization of power generation in the driver 
fuel (because all such heat must be rejected to a sink) by varying dimensions and mate- 
rials. A reference configuration was obtained that consisted of a spherical test cavity 
surrounded by driver fuel elements arranged in a cylindrical annulus. A heavy water 
moderator region between the cavity and the driver fuel provided thermalization of the 
neutron source to the cavity. The reactor was reflected both radially and axially with 
25. 4 centimeters of beryllium. Driver fuel loading was 6. 4 kilograms of enriched 
uranium, and driver power generation was 41. 3 megawatts. Required control swing 
for 1200 -megawatt hours of operation was 18 percent reactivity Ak/k. This could be 
adequately provided by steel-clad cadmium hollow control rods. Twelve rods containing 
0. 152-centimeter-thick cadmium provided 50 to 60 percent Ak/k when interspersed in 
the driver fuel region. Ancillary results of material reactivity worths, reactivity coef- 
ficients, flux spectra, and power distributions were reported. 


INTRODUCTION 

Experimental and analytical work on a coaxial-flow, open-cycle gas-core nuclear 
rocket concept has been in progress for 10 to 15 years (refs. 1 to 3). This concept is 
basically a critical mass of nuclear fuel in the form of a plasma ball suspended by fluid 
dynamic forces in the center of a flowing propellant. The propellant is contained in a 


reflector and pressure vessel assembly. Since all structural parts are thermally 
shielded by the propellant (via an injected seed material), fuel temperatures on the 
order of 50 000 K can be theoretically achieved. Consequently, propellant temperatures 
required to obtain specific impulses around 2000 to 5000 seconds are theoretically pos- 
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sible. The concept also has the capability of producing thrusts on the order of 10 to 
10° newtons. This combination of high specific impulse and high thrust would represent 
a major breakthrough in rocket propulsion technology if sufficiently low engine weights 
and fuel loss rates could be achieved. 

As part of recent gas-core studies, a cavity test reactor has been proposed (ref. 4). 
This facility would be used to experimentally test the fundamental phenomena upon which 
the gas-core rocket concept is based. Conceptually, the test reactor would be a con- 
ventional solid-core reactor with a large test chamber in the center. Thus, a driver 
core would be used to provide neutron flux to a plasma fuel in the test chamber. In this 
manner the plasma fuel criticality requirement is removed, thereby significantly reduc- 
ing core size and pressure that would be required for a gas-core reactor. 

This report presents the neutronics information generated from a preliminary de- 
sign study of such a cavity test reactor. Although a feasible configuration was obtained 
in the study, the final design was not optimized. The design presented herein was aimed 
primarily at minimizing the reactor power level within the constraints of the study. An 
optimized design would require a tradeoff between operating and capital costs, and the 
necessary information for such a determination was not obtained during the study. 

Major constraints of the study were that existing fuel-element technology be used in 
the driver core and that the capability exist for tests in the cavity. Thus, the driver 
core was designed around MTR fuel elements, thereby using considerable reactor oper- 
ating experience at the Plum Brook Reactor. Also, the cavity will be 60. 96 centimeters 
(2 ft) in diameter. Under test conditions (calculated from methods described in ref. 5) 
the cavity will contain 375 grams of enriched uranium, which will generate 2. 7 mega- 
watts of power. All configurations will be restricted to single fluid systems; that is, if 
a liquid moderator proves desirable, that material will also be used to cool the fuel ele- 
ments and the reflector. 


ANALYSIS 

Objective 

The neutronics design effort was directed at varying materials and configurations in 
order to minimize the reactor power generated in the driver fuel region while maintain- 
ing the required power in the test cavity. All driver power must be rejected to a heat 
sink and will be subsequently referred to as throwaway power. 
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Since the cavity power was fixed at 2. 7 megawatts, the power split Pq/Pq or the 
ratio of driver power to cavity power, can be considered a measure of throwaway 
power. The power split can be represented by 
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where 

Mp driver fuel mass, kg 

M^, cavity fuel mass equal to 0. 375 kg of uranium 235 
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Of(E) energy dependent microscopic fission cross section, cm 

— 2 
cp(E, V) energy and space dependent neutron flux, neutrons/cm -sec 

3 

driver fuel volume, cm 

3 

cavity fuel volume, cm 

To simplify equation (1) we can assume that the flux spectra are the same in both fuel 
regions and that all power is produced by thermal fissions. Then 


P D _ M P^th, D ( 2 ' 

P C 0-3T5cP th)C 

Equation (2) indicates that in order to minimize Pj^/Pq, M d should be minimized 
and/or the ratio of thermal fluxes in the driver and the cavity regions reduced. Another 
variable to consider is the fast flux from the driver fuel region c p, which is the 
source for £,. 

Unfortunately, all these quantities are interdependent, and physical changes in the 
reactor design tend to have offsetting effects on P^/P^. For example, a moderating 
reflector tends to lower but reduces cp ^ which in turn increases cp^ p/^th C‘ 
However, there does exist some combination of reactor materials and configuration that 
will certainly tend to minimize Pj)/Pq- In general one can expect that hardening the 
spectrum in the driver region will increase and lower jyVj-fo q- Also, for a 

given moderator material some optimum thickness exists that will balance the slowing 
down and absorption effects on <p^ -q- 
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Calculational Model 


To satisfy both mechanical design and experimental test requirements, the reactor 
consists of a spherical test cavity surrounded by a neutron moderating region, a cylin- 
drical driver fuel region, reflector, and pressure vessel. A description of the reactor 
and local cavity conditions are shown in figure 1. The cavity conditions are those re- 
quired for a gas-core reactor test. The thickness of the driver fuel region is fixed by 
the selection of an MTR fuel element used in the Plum Brook reactor (fig. 2), although 
the absolute location of this region is variable. Top and bottom reflectors were con- 
sidered necessary both to reduce driver fuel mass and to produce an even flux distribu- 
tion in the test cavity. 

For initial parametric calculations a spherical model was used in which the various 
components were homogenized and arranged in spherical shell regions (fig. 3). How- 
ever, because a spherical shell is a considerable deviation from the cylindrical arrange- 
ment of driver fuel elements, more accurate values for multiplication factors and flux 
levels were obtained from two-dimensional models in RZ (cylinder) and R0 (disk at 
reactor midplane) geometry. The RZ model (fig. 4) required that the driver fuel re- 
gion be homogenized into a cylindrical annulus and that the cavity regions be cylindrical. 
In R0 geometry (fig, 5) it was possible to explicitly describe the driver fuel elements, 
but axial dimension effects had to be simulated by the use of axial leakages generated in 
an RZ calculation. 

As noted previously, each calculational model contained certain restrictions. The 
best values for multiplication factor were considered to be those from the R0 calcula- 
tion with appropriate correction factors for energy groups and cross section detail. With 
the use of zone dependent axial leakages, the Rd model could account for three- 
dimensional heterogeneity effects in the driver fuel elements. 

Whenever the same reactor configuration was represented by different calculational 
models, dimensions were adjusted to conserve material volumes. This was considered 
a better technique than attempting to preserve mean free paths across regions. 

The driver fuel height was set at 86. 6 centimeters (cylindrical model) to provide a 
reasonable viewing angle from the test cavity fuel region. Hopefully, this geometry 
would provide an equal source of neutrons to the entire test cavity fuel region. No at- 
tempt was made to optimize the fuel height during this study. 


Reactor Codes 

All reactor calculations were performed with multigroup neutron transport codes; 
TDSN (ref. 6) for one-dimensional calculations and DOT (ref. 7) for two-dimensional 
calculations. Generally, the one-dimensional (spherical) calculations were in S^P^IQ 
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group detail, and the two-dimensional (cylindrical and R0) calculations were in SgPglO 
group detail. Correction factors were calculated so that two-dimensional S^P^IS 
group results could be synthesized. 

Cross sections were generated with GAM-n (ref. 8) and GATHER-II (ref. 9) codes 
for energy groups above and below 2. 38 electron volts, respectively. For the purpose 
of cross-section generation all groups above 2. 38 electron volts were treated as fast 
neutrons and those below 2. 38 electron volts were treated as thermal neutrons. Cross 
sections were flux weighted for the particular material regions in which they were to be 
used. The 19 -group energy set (12 fast and 7 thermal) had been developed and used in 
previous analytical and experimental work (ref. 3). Full down- scattering was allowed 
between the fast groups and full down- and up-scattering were allowed between the ther- 
mal groups. The 10-group energy set (five fast and five thermal) was merely a con- 
solidation of the 19-group structure. The same scattering provisions were retained in 
the 10-group energy set. In the high temperature regions in the test cavity free atom 
scattering kernels were used to determine temperature effects on cross sections. 


Calculation Interpretation 

For the most part the calculated configurations were noncritical and contained two 
fuel zones, which were neutronically coupled. Thus, each calculation produced a mul- 
tiplication factor for the test cavity k^ and one for the driver fuel k D ; the sum of 
which represented the multiplication factor of the reactor k (see appendix A). Reac- 
tivity effects for noncritical systems were determined from 


Ap = p 2 - Pl 


k 2 ~ k l s Ak 
k l k 2 k 


(3) 


where p is reactivity and the subscripts 1 and 2 are the initial and perturbed states, 
respectively, of the reactor. 

For a reactor with two fuel zones, local reactivity effects were calculated from 


Ak C s k C2 ~ k Cl 
k k i k 2 


and 



Ak D _ k D2 " k Dl 


( 5 ) 


k k l k 2 


The use of reactivity effects to adjust the multiplication constants of a calculated con- 
figuration is facilitated by manipulating equations (4) and (5) to obtain 



and 



Equations (6) and (7) can then be used to calculate k D and of a perturbed reactor 
if the initial conditions and the reactivity effect of the perturbation are known. 


RESULTS AND DISCUSSION 
Fuel Loading 

Neutronic coupling of the two fuel zones can be observed by varying the fuel loading 
of the driver fuel region M Q . Data plotted in figure 6 show k D increasing and k^ 
decreasing with increasing Mp. Increased neutron absorption caused by the addition of 
driver fuel tends to decrease cavity flux thereby decreasing k^,. Since P^/Pq ^kp/k^, 
the conclusion is obvious that for a given configuration the best power split is obtained 
with the lowest driver fuel loading. 

Specific fuel reactivity worths can be derived from figure 6 by 

F= - Ak /L (8) 

AM/M 
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where AM/M represents the change in fuel mass divided by the average fuel mass. 
Values for the test cavity and driver fuel regions are presented in table I and figure 7. 
These data are useful in adjusting the fuel loading of a calculated configuration to obtain 
a given multiplication constant. 


Moderator Thickness 

Materials considered for the moderator were HgO, beryllium, D 2 O (0. 25 vol. % 
H 2 O). For each moderator material some optimum thickness exists for which throw- 
away power is a minimum at a given cavity power. Thinner moderators would not pro- 
vide sufficient neutron thermalization, and thicker regions would cause excessive neu- 
tron absorption. In either case fewer fissions per unit source occur in the cavity test 
region. Also, changes in driver fuel volume resulting from varying moderator thickness 
produce a second-order effect on power split. At smaller fuel volumes (smaller mod- 
erator thickness) driver fuel mass tends to decrease thereby decreasing the power split. 
This occurs because the driver fuel region approximates an infinite slab with constant 
thickness; therefore, its fuel density required for criticality is nearly constant. 

Calculated results show throwaway power to be quite sensitive to H 2 0 thickness 
(fig. 8). A minimum of 65 megawatts occurred at 2. 54 centimeters (1 in. ). Values of 
99 megawatts were obtained at the lower limit (no moderator) and 200 megawatts at 
7. 62 centimeters (3 in. ). The ratio for DgO was lower than for HgO and was 

relatively flat as a function of moderator thickness, although a minimum of 52 megawatts 
did occur at 15. 24 centimeters (6 in. ). One point at a 15. 24-centimeter (6-in. ) thick- 
ness was calculated for beryllium (DgO cooled), which indicated a Pj-j/Pq of 60 mega- 
watts, somewhat higher than for DgO. In all cases the driver fuel was cooled with the 
moderator material (except when Be was used). Greater absorption in the HgO fuel 
coolant contributes to the HgO minimum being greater than the DgO minimum. 

Heavy water was selected as the moderator material for two reasons: its lower 
throwaway power (and therefore P^/P^) and its insensitivity to thickness, which was 
considered desirable from a mechanical design standpoint because the moderator-region 
has a spherical inner surface and a cylindrical outer surface. This selection also fixed 
D 2 O as the coolant for the driver fuel and the reflector because of the decision to main- 
tain a single fluid system. 

The effect of the moderator on reactivity is assumed to be directly related to total 
collisions. Therefore, in a highly moderated reactor where the scattering is nearly 
isotropic, this effect will be retained between calculational models by keeping the mod- 
erator volume constant when relating a spherical model to the reference configuration. 
Thus, the minimum thickness of 15 centimeters from figure 8 was not as significant as 
the corresponding volume of 3.9x10^ cubic centimeters. 
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Since the purpose of figure 8 is to show the relative effect of moderator material 
and thickness, data based on an early model of the cavity test reactor were plotted. 

These data are not consistent with the final configuration because the cavity contained 
559 grams of uranium-235 at 2200 K and hydrogen propellant at 5000 K with 17 weight 
percent uranium-238 seed. Also, the results are from spherical calculations and have 
not been corrected to a cylindrical model. 

The magnitude of relative effects from the use of different moderator materials was 
observed to be sensitive to the particular configuration. For example, if coolant-filled 
control rod channels were included in the driver fuel zone, the resulting change in region 
composition would affect the throwaway power comparison of the HgO and D 2 O moder- 
ated reactors. Since a greater volume fraction of the fuel zone would be occupied by 
coolant, the effect of changing coolants would be magnified. 


Reflector Thickness 

Materials considered for the reflector were beryllium and DgO. The decision to 
use D 2 O as moderator eliminated HgO as a possibility because of the single-fluid 
groundrule. The primary effect of the reflector appears to be a reduction of driver fuel 
critical mass. A secondary effect is to provide a portion of the neutron source to the 
test cavity by reflecting fast neutrons that would penetrate the driver fuel region. As 
reflector thickness increases, critical mass decreases, which tends to decrease throw- 
away power. The primary effect is enhanced by a moderating reflector material 
whereas the secondary one requires a fast neutron reflector (generally a high molecular 
weight material). This would indicate that some material with intermediate properties 
and/or a composite reflector would be desirable. In addition, absorption of neutrons in 
the reflector should be minimized to decrease driver fuel mass required for criticality. 

Calculated results are presented in figure 9 for various reflector thicknesses, all 
of which are backed by a 10. 16 -centimeter (4-in. ) steel pressure vessel. Reflector 
thickness refers to a cylindrical radial reflector plus top and bottom slabs, all of equal 
thicknesses. Spherical model reflector thicknesses have been converted from this model 
by preserving total reflector volume. All results have been normalized to a cylindrical 
model of the reference configuration with a test cavity power of 2. 7 megawatts. The 
lowest throwaway power Qp, 40. 1 megawatts, was obtained with the largest beryllium 
reflector thickness (35. 6 cm (14 in. )) that was calculated. The throwaway power 
increased at decreasing beryllium thickness to a value of 49. 4 megawatts at 10. 16 cen- 
timeters (4 in. ). The throwaway power for D 2 O reflectors was about 16 megawatts 
higher for equivalent thicknesses. One composite reflector (8. 89 cm Be + 12. 45 cm 
D 2 O) resulted in a 3 -megawatt increase in Qp for the same total thickness. 
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Based on these results, a 25. 4-centimeter (10-in. ) thick beryllium reflector was 
selected for the reference configuration. On an equivalent thickness basis beryllium 
was clearly the better material. However, the particular thickness selected was some- 
what arbitrary. For a final selection, cost data would be required to evaluate the trade- 
off between capital cost of the material, and the operating cost of rejecting Q^. Also, 
if the beryllium material costs are particularly high, D2O could be used to supplement 
a thinner beryllium reflector, or perhaps even eliminate it. 

The fast neutron reflectivity contribution of the pressure vessel was not explicitly 
determined in these calculations. Any change in pressure vessel thickness or material 
might have an effect on the results of the thin beryllium reflectors and on the DgO re- 
flectors. It is doubtful, though, that any effect would be noticeable for the case of the 
relatively thick beryllium reflected configurations considered. 


Reference Core Characteristics 
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Based on data presented in this report, a DgO moderator with a volume of 3. 9x10 
cubic centimeters and a 25. 4-centimeter-thick beryllium reflector was selected. These 
items were included with the nonnuclear test dependent design features to obtain the 
reference model configuration of the cavity test reactor (table n). The neutronic charac- 
teristics of this reference reactor (itemized in table m) indicate that a driver fuel mass 
of 6. 4 kilograms uranium (0. 932 U) is required to provide excess reactivity for 
1200 megawatt -hours of operation. Other data include median fission energies of 
0. 06 and 0. 04 electron volt in the test cavity and driver fuel regions, respectively. The 
higher value in the test cavity can be attributed to upscattering of neutrons in the high 

temperature propellant and fuel plasma regions. A power level of 41.3 megawatts in 
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the driver region is required to produce the required flux level of 4. 1x10 neutrons per 

square centimeter per second to the cavity to generate 2. 7 megawatts. 

Fast (E > 0. 07 MeV) and thermal (E < 0. 12 eV) flux distributions along radial and 

axial midplanes of the reactor are plotted in figure 10. Of particular note is the flat 

distribution of thermal flux in the test cavity which predicts a flat power distribution. 

As expected the thermal flux showed peaks in the moderator materials and depressions 

in the absorbers. The fast flux peaked in the fuel regions and tended to drop off rapidly 

in the axial reflector. Although not shown in figure 10(b), the reduction from the center 

-4 

to the upper edge of the beryllium reflector was about 10 . 

Median fission energies of 0. 04 and 0. 06 electron volt calculated for the driver and 
cavity fuel regions indicate that the reference model is a thermal reactor. Flux spectra 
for the two fuel regions are plotted in figure 11. Displacement of the slow neutron en- 
ergy distribution in the cavity toward higher energy (compared with the driver fuel re- 
gion) substantiates the calculated higher median fission energy for that region. Moder- 
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ation in the D 2 O tends to decrease the cavity spectrum at energies below the fission en- 
ergy peaking (~1. 6 MeV) until the upscattering effect becomes apparent around 2 elec- 
tron volts. 

Distribution of power generation in the test cavity fuel is important to heating rates 
and associated flow effects. Peak-to-minimum values of 1. 1 for the power density are 
indicative of a flat distribution (fig. 12). Power density in the axial direction is slightly 
greater than along the radial midplane even though no fuel is located above or below the 
cavity. Apparently the axial reflector region acts as a flux trap for thermal neutrons 
and provides a source for the cavity region. 

Variations from the average in figure 12 are somewhat deceptive because the linear 
dimension represents a radius and the power density is volume averaged. Also, it 
should be noted that these data are based on a homogeneous distribution of fuel atoms 
because of a lack of any experimentally predicted distribution for reference model con- 
ditions. 

Power density traverses through the driver fuel elements in all three dimensions 
are presented in figure 13. The hot channel in a fuel element occurs along the inner 
plate (moderator side) where the local-to -average power ratio is about 1. 3 (fig. 13(a)) 
at the axial average position. Power also peaked on the side of a fuel element adjacent 
to an empty control rod channel (occupied by DgO). The maximum hot spot on that sur- 
face was 1. 8 (fig. 13(b)) at an average axial position. Axial power distributions have 
been averaged and presented as relative values in figure 13(c). Therefore, multiplica- 
tion by local-to-average ratios in figure 13(c) can be used to convert power density val- 
ues in figure 13(a) and (b) to a specific axial location. The expected axial power shape 
(cosine distribution) was obtained with a sharp upturn near the edge caused by thermal 
neutrons from the beryllium reflector. 


Control System 

Excess reactivity requirements are based on reactor startup under nontest condi- 
tions and with the test cavity flooded with hydrogen at 300 K and 200 atmospheres. The 
negative reactivity effect of that hydrogen (-5. 7 percent Ak/k) was the major component 
in the total 9 percent Ak/k required (table IV). Other items were fuel depletion, 

0. 3 percent Ak/k; °Xe production, 0. 4 percent Ak/k; other fission products, 1 per- 
cent Ak/k; D 2 O temperature defect, 0. 1 percent Ak/k; and an added contingency of 

1. 5 percent Ak/k. Fuel depletion was based on total absorption of neutrons by fuel 
atoms during the 1200 -megawatt-hour operation. The xenon- 135 penalty was calculated 
for a 30-minute test run. The temperature defect results from the average D 2 O tem- 
perature rising from 300 to 325 K at full power. A contingency was added to account for 
calculational and design uncertainties. 
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Although the worth of 375 grams of uranium (0. 932 U) is about 7 percent Ak/k 
(or about 10 dollars) a corresponding negative effect of -3 percent Ak/k occurs in the 
driver fuel. Therefore, the net effect of adding fuel to the cavity for a test run is 4 per- 
cent Ak/k. The control system must have sufficient capacity to override this effect and 
to shutdown the reactor under test conditions with an estimated safety margin of 5 per- 
cent Ak/k. Consequently, total required control swing is 18 percent Ak/k, and the 
corresponding k's for excess reactivity and shutdown are 1.099 and 0.917, respec- 
tively. 

Three control system configurations were proposed to obtain this control swing. 

All reference configuration calculations were based on MTR control rods (type A) in- 
serted directly into every third driver fuel element position, for a total of 12 rods in the 
reference configuration. These rods were hollow, 7.62-centimeter-square rectangular 
boxes. The walls were 0. 952-centimeter-thick stainless steel with a 0. 1524-centimeter- 
thick sheet of cadmium (Cd) sandwiched in the stainless steel. The hollow portion of the 
rods was occupied by a coolant. Types B and C were hollow tubes with 6. 9 8 -centimeter 
outside diameters. Materials and their thicknesses were the same as type A. These 
tubes were inserted into a stationary channel formed by a zircalloy-4 tube with a 0. 635- 
centimeter-thick wall and a 7. 62-centimeter outside diameter. Types B and C differed 
only in the control rod coolant, I^O for type B and a gas for type C. 

Control swings were obtained from spherical calculations using smeared control 
rod regions and a spatial self- shielding factor of 0. 02 for the cadmium (appendix B). 
Reactivity differences for rods in and rods out calculations indicated a swing of about 
50 percent Ak/k could be obtained for any of these configurations (table V). In order to 
increase confidence in the spherical calculations, a two-dimensional calculation in R0 
geometry was performed in which control rods and fuel elements were explicitly defined. 
Control rod cross sections were obtained from a one -dimensional cell calculation. The 
resulting control swing was 66 percent Ak/k, thereby indicating that control swing pre- 
dictions from the spherical calculations are probably conservative. Thus, either type 
A, B, or C would provide adequate control swing. In a more detailed design excess 
swing could be reduced by varying the cadmium content. Also, the rods would be as- 
signed various functions, such as scram, shim, and regulating, and the poison content 
and/or drive mechanism would be adjusted accordingly. 

In addition to control swing, the control system configuration also affects the mul- 
tiplication factor of a reactor. For comparison purposes variations in k are translated 
to required power generation in the driver fuel by normalizing all systems to a test 
cavity power of 2. 7 megawatts. Thus, the reference power of 41. 3 megawatts is in- 
creased to 48. 1 and 42. 9 megawatts when type B and C control systems are used (ta- 
ble VI). Added neutron absorption from inclusion of the zircalloy-4 tubes and the HgO 
coolant could account for this increase. 
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Reactivity Coefficients 


For use in reactor dynamics analyses, a number of reactivity effects were obtained 
from static calculations of various design perturbations. The data are presented in ta- 
ble VII explicitly for each fuel zone of a reactor under test conditions. Reactivity 
changes in each fuel region were calculated from equations (4) and (5). In general, op- 
posing reactivity changes occurred in the two fuel regions when an operating condition 
was perturbed. Such changes should be amenable to regulation by the control system. 
Insertion of control rods into the driver fuel region tended to reduce reactivity in both 
fuel zones at a relative rate k^/k^, of about 100 to 1. 

These results represent a partially self -regulating system from the standpoint of 
overall reactor control. However, an additional function of the control system is the 
maintenance of sufficient flux in the test cavity to produce 2. 7 megawatts. 

Perturbations used in table VII are relatively straightforward except for propellant 
temperature, in which only the effect of temperature on microscopic cross section was 
considered. Hydrogen atom density was not changed corresponding to the temperature 
change. The effect of hydrogen atom variation can be derived from the calculated pro- 
pellant pressure coefficient. 

Material reactivity worths determined during the course of the study are itemized 
in table VIII. Of particular note is the small penalty incurred if tungsten (W) is used 
as the propellant seed material (since W is a high absorbing material) and/or if the 
D 2 0 is contaminated to 1 percent of its volume with HgO. To maintain the 2. 7-megawatt 
test cavity power, the use of W would require a 0. 8-megawatt increase in throwaway 
power and contamination by 1 percent HgO, a 1. 3-megawatt increase. Addition of hy- 
drogen uniformly to the fuel plasma in the test cavity significantly reduced local reac- 
tivity thereby requiring increased throwaway power generation. If the 375 grams of 
cavity fuel were to deposit (at the plasma density) on the inner surface of the cavity 
liner, cavity reactivity would increase 2 percent Ak^/k but a negative effect of 1 per- 
cent Akjj/k in the driver fuel tempers the overall effect. 

SUMMARY OF RESULTS 

A neutronic analysis of a cavity test reactor, designed to test the feasibility of the 
gas-core reactor concept, was performed. The spherical cavity test section is 60. 96 
centimeters in diameter and is surrounded by a neutron moderator region, a cylindrical 
driver fuel region containing MTR fuel elements, a neutron reflector, and a pressure 
vessel. Primary consideration in component material selection and sizing was given 
to minimizing throwaway power; that is, power generated in the driver fuel in order to 
provide sufficient flux to the test cavity to produce the 2. 7 megawatts necessary for gas- 
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core reactor experiments. Thus, heavy water (D 2 O) was selected to moderate the 

source flux to the cavity and to cool the driver fuel elements. Moderator volume is 
5 

3. 9x10 cubic centimeters. The reflector is a 25. 4-centimeter-thick cylindrical annulus 
of DgO cooled beryllium with 25. 4-centimeter-thick end pieces. A 10-centimeter-thick 
steel pressure vessel is required to contain the operating pressure of 200 atmospheres. 
Twelve control rods and 24 fuel elements are arranged in a 8. 4-centimeter-thick cylin- 
drical annulus with an inner radius of 45. 94 centimeters. This driver fuel region con- 

235 

tains 6. 4 kilograms of enriched uranium (0. 932 U) and generates 41. 3 megawatts. 
Specific results of the study are the following: 

1. Throwaway power was increased by 25 and 14 percent (compared with DgO), when 
water (HgO) and beryllium, respectively, were used as moderator materials, primarily 
because of greater parasitic absorption of neutrons. These numerical values are sensi- 
tive to the assumed configurations, in particular, the volume fraction of coolant in the 
driver fuel region. 

2. For the reflector thicknesses of interest (<36 cm), a lower throwaway power is 
obtained using a beryllium reflector than a D 2 O reflector because of reduced critical 
mass and reflectivity of fast neutrons. These neutrons penetrate the driver fuel region 
and contribute to the flux in the test cavity. 

3. An excess multiplication factor of 1. 099 for a clean reactor with a void test cav- 
ity is required based on 1200-megawatt-hour operation with an outlet coolant tempera- 
ture of 325 K. Total required control swing is 18 percent Ak/k, which includes a 5 per- 
cent Ak/k shutdown margin. 

4. Fuel required for a gas-core reactor test (375 g U) is worth 7 percent Ak/k 
locally but only 4 percent Ak/k for the overall reactor. This results from neutronic 
coupling between test cavity and the reactor fuel zones. 

5. Hollow rod control rods containing cadmium can provide a control swing of about 
50 to 60 percent Ak/k. A control rod is located in every third position in the 

36 -element annular fuel channel. 

6. Reactivity coefficients were calculated for various system perturbations. In 
general, a system change affects the reactivity of the cavity test fuel and the reactor 
driver fuel in opposite directions. 

7. Required power output was insensitive to the selection of seed material and to 

contamination of D„0 with H 0 0 (up to 1 percent). 

^ ^ 14 

8. Average thermal flux (E < 0. 12 eV) in a void test cavity is 4. 1x10 neutrons 

14 

per square centimeter per second and at gas -core test conditions is 1. 6x10 neutrons 
per square centimeter per second. Neutron spectra are near thermal as evidenced by 
median fission energies of 0. 06 and 0. 04 electron volt in the test cavity and driver 
reactor fuel regions, respectively. 
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9. Power distribution in the test cavity is relatively flat with a peak- to -minimum 
value of 1. 1. In the reactor fuel elements, three-dimensional power distributions were 
calculated which showed local peaking at the axial center, in the inner fuel plate, and 
along the side next to a coolant filled control rod channel. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 22, 1972, 

503-04. 
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APPENDIX A 


CALCULATED MULTIPLICATION CONSTANTS 

Comparison of the various design configurations was made on the basis of required 
power output in the driver fuel to obtain 2. 7 megawatts of power in the test cavity. To 
do this, all designs were normalized to a k D that would provide the required reactivity 
margin. 

Based on fuel depletion, temperature defect, cavity hydrogen flooding, and fission 
product buildup, an excess reactivity of 9 percent Ak/k is required for a reactor under 
nontest conditions (no fuel or hydrogen in the test cavity). The corresponding k is 
1.099. Under the assumption that an R0 S^Pjl9 group result using RZ calculated 
axial leakages represents the best accuracy that is obtainable within the practical limits 
of computer capacity and running time, correction factors are used to adjust the k from 
the spherical calculations used in the parametric analyses. To reduce the number of 
calculations needed, the related k^ for a system under test conditions is used for nor- 
malization. In effect, normalizing the results to this k^ should provide the best cal- 
culational accuracy (within the limits of the study) and a basis for comparing different 
design configurations. 

Calculated values for adjusting multiplication factors are listed in table IX. Inclu- 
sion of fuel and hydrogen in the test cavity significantly reduced k D (by 3. 18 percent 
Akp/k) because of increased absorption of neutrons, which had the potential of causing 
fissions in the driver fuel. This is another illustration of neutronic coupling between 
the fuel zones. Because of this coupling, adjustments of k D must be made with equa- 
tion (6). Since cylindrical calculations were limited by machine capacity and running 
time to S 2 Pq 10 group models without pressure vessels, a correction factor was calcu- 
lated using spherical models. The R0 -to- spherical model correction of 5. 48 percent 
Akp/k was primarily due to changing the driver fuel from a cylindrical arrangement to 
a spherical shell. Application of all correction factors leads to the result that a spher- 
ical S^Pq 19 group calculation of k D “ L 135 for a reactor under gas core test condi- 
tions would have a "best calculation" value of 1. 099 for a reactor with a void cavity. 

Variation between calculational models also occurs in k^. This is accounted for 
by comparing k^/k^, for two-dimensional and spherical calculations. Thus, spherical 
power splits Pj)/Pq = k^/k^, are multiplied by 0. 956 to obtain the more accurate two- 
dimensional results. 

In practice a ID S^P^19 group calculation is performed on a given configuration; 
kj-j is normalized to 1. 135 by adjusting the fuel loading (fig. 7). The term k^, is ad- 
justed, because of the change in fuel loading, using figure 7 and equation (7). Then the 
ratio of k^/k^, is multiplied by 0. 956 and driver power Pj-j is calculated based on 
2. 7 megawatts cavity power Pq. 
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APPENDIX B 


SPATIAL SELF SHIELDING IN CADMIUM 

For use in smeared control rod regions used in spherical calculational models, a 
self shielding factor f was calculated in order to account for the thermal flux depres- 
sion in the cadmium section of the control rod. In theory, f is needed because a similar 
flux depression would not exist in a smeared region. Therefore, the cadmium cross 
section must be adjusted by f so that the total reaction rate will be preserved. 

It was assumed that the rectangular channel geometry of an MTR control rod could 
be approximated by a slab model. For large values f = l/2Et where S is the total 
cross section in barns and t is the thickness in centimeters of the cadmium (ref. 10). 
Energy dependent values for f were determined to be about 0. 02 for the energy range 
where most reactions are expected to occur (table X). Therefore, a single value of 
0. 02 was used for all energy levels in the study presented herein. Based on later cal- 
culations with cell averaged cross sections, a somewhat higher value for f would have 
been more accurate. 
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TABLE I. - SPECIFIC REACTIVITY WORTH 


OF DRIVER FUEL 


Average 

driver 

fuel, 

hg 

Mass 

ratio, 

AM/M, 

percent 

Cavity re- 
activity 
effect, 
Ak c /k, 
percent 

Driver re- 
activity 
effect, 
Ak D /k, 
percent 

Specific fuel reac- 
tivity worth a , F 

In cavity 

In driver 

5. 85 

15. 39 

-0. 196 

2.97 

-0.0126 

0. 193 

6. 75 

13. 33 

-. 146 

2. 15 

-. 0110 

. 161 

7. 65 

11. 76 

120 

1.75 

-.0102 

. 149 

8. 56 

10.75 

-.091 

1. 42 

-.0085 

. 132 


a F = (Ak/k)/(AM/M); where M = (M, + M 2 )/2. 


TABLE II. - REFERENCE MODEL CONFIGURATION 


Component 

Material 

Shape 

Radius, 

cm 

Height, 

cm 

Wall thick- 
ness, 
cm 

Test cavity liner 
Hydrogen plenum 

Aluminum 

Sphere 

Sphere 

30. 48 

31. 115 


0. 635 
2. 54 

Moderator inner liner 
Moderator 

Aluminum 
Heavy water (DgO) 

Sphere 

33. 655 
a 34. 29 
b 45. 315 

34. 29 
43. 3 

.635 

Moderator outer liner 

Zircalloy-4 

Cylinder 

45. 315 

43.3 

.635 

Driver fuel channel 

MTR fuel elements 
and DgO coolant 

Cylinder 

a 45. 94 
b 54. 347 

43.3 


Fuel channel liner 

Zircalloy-4 

Cylinder 

54. 347 


. 635 

Reflector 
Pressure vessel 

Beryllium and DgO 
coolant 
Steel 

Cylinder 

Cylinder 

a 54. 982 
b 80. 382 
80. 382 

a 43. 935 
b 69. 335 

10. 16 


a Inner. 

b Outer. 





TABLE III. - NEUTRONIC CHARACTERISTICS OF 


CAVITY TEST REACTOR 


Test cavity fuel loading a , kg U 

0. 375 

Driver fuel loading a , kg U 

6.4 

Average thermal neutron flux in test cavity* 3 , neutrons/cm^ sec 

4. lxlO 14 

For nontest conditions (void cavity) 

For operating test conditions 

1.6X10 14 

Median fission energy, eV 


Test cavity fuel 

0. 06 

Driver fuel 

0. 04 

Power generation, MW: 


Test cavity fuel 

2.7 

Driver fuel 

41. 3 


a Atomlc fraction of u, 0. 932. 

b E ^ 0. 12 eV. 


TABLE IV. - REACTIVITY REQUIREMENTS FOR 


1200 MILLIWATT-HOUR OPERATION 



Reactivity, 

Ak/k 

Excess: 


Fuel depletion 

0.003 

Xenon- 135 production 

.004 

Other fission production (estimated) 

.010 

Temperature defect (300 to 325 K DgO) 

.001 

Hydrogen flooding of test cavity 

.057 

Subtotal 

. 075 

Contingency 

. 015 

Total 

.090 

Shutdown: 


Addition of test cavity fuel 

.04 

Margin allowance 

. 05 

Total 

. 09 

Required control swing 

. 18 

Effective multiplication factor (void test cavity), k^ 

1.099 

Shutdown multiplication factor, k g 

. 917 
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TABLE V. - CALCULATED CONTROL SWING FOR VARIOUS 


CONTROL ROD CONFIGURATIONS 


Type 

Description 

Coolant 

Cadmium 

surface 

area, 

cm 2 

Reactivity 1 , 

Ak/k 

A 

MTR type, 7. 62 cm 2 

Heavy water 

4. 6xl0 3 

0. 49 
b . 66 

B 

Hollow tube, 6. 98-cm o. d. 

Water 

3. 3xl0 3 

. 57 

C 

Hollow tube, 6. 98-cm o. d. 

Gas 

3. 3xl0 3 

. 52 


a Spherical calculations with smeared control rod regions and a cadmium 
spatial self-shielding factor of 0.02. 

^Two-dimensional calculation using cell averaged cross sections for the 
control rod region. 


TABLE VI. - REQUIRED POWER AS FUNCTION OF 
CONTROL SYSTEM CONFIGURATION 


Control 

Description 

Driver fuel 

system 


power, 



MW 

A 

MTR type, heavy water cooled 

41. 3 

B 

Tube type, water cooled 

48. 1 

C 

Tube type, gas cooled 

42. 9 





TABLE VII. - REACTIVITY EFFECTS OF DESIGN PERTURBATIONS 


IN THE REFERENCE CONFIGURATION 


Variable 

Range 

Specific fuel re- 


Lower 

Upper 

activity worth 


F c 

f d 

Test cavity fuel, kg 

0. 3375 

0. 375 

0. 035 

-0.017 

Test cavity fuel expansion, cm 

19. 2 

20. 42 

. 023 

-.011 

Propellant pressure, A 

160 

200 

017 

. 0094 

Propellant temperature, K 

2560 

5000 

-. 0186 

. 0183 

Plenum pressure, A 

160 

200 

-. 0015 

-. 0118 

DgO pressure, A 

1 

200 

0 

3. 5x10'® 

D 2 O temperature, K 

300 

325 

0 

-. 0097 


a F = (Ak/k)/Ap/p where Ap represents the change in the particular 
property and p is the arithmetic average of the lower and upper 
limits. 


TABLE VIII. - MATERIAL REACTIVITY WORTHS IN 
THE REFERENCE CONFIGURATION 


Variable 

Cavity re- 
activity 
effect, 
Ak^/k 

Driver re- 
activity 
effect 
Akjj/k 

Change in 
driver 
throwaway 
power, 

AQd, 

MW 

Addition of 12.6 g of hydrogen uniformly to 

-0. 0121 

0. 0069 

8.7 

test cavity fuel 




Redistribution of test cavity fuel onto inner 

.020 

-. 010 

--- 

surface of cavity liner 




Substitution of W in place of ^®®U seed 

-. 0004 

-.001 

.8 

Inclusion of 1% HgO in DgO 

-. 0005 

-. 0061 

1. 3 





TABLE IX. - REACTIVITY CORRECTIONS FOR ADJUSTING CALCULATED 


MULTIPLICATION CONSTANTS 


Variable 

Driver re- 
activity 
effect, 
Ak D /k 
percent 

Cavity re- 
activity 
effect, 
Akj-./k 
percent 

Driver multi- 
plication 
factor, 

k D 

Initial reactor with no cavity fuel or 


— 

1. 099 

hydrogen 




Inclusion of 375 g uranium and hydrogen 

-3. 18 

4.70 

1. 060 

in cavity at test conditions 




SgPglO group (no pressure vessel) to 

. 24 

-. 02 


S^P^19 group model 




RS model a to spherical model 

5. 48 

. 10 

1. 135 


3 . 

Using axial leakages from an RZ cylindrical calculation. 


TABLE X. - SELF-SHIELDING FACTOR 
FOR 0. 1524-CENTIMETER-THICK 


CADMIUM SLAB 


Upper energy 
of group, 
eV 

Total cross 
section, 
b 

Self- shielding 
factor, 
f 

2. 38 

21. 1 

0. 16 

. 414 

56. 2 

. 058 

. 12 

168. 6 

. 0281 

.08 

117. 0 

. 0281 

. 0253 

137. 1 

. 0239 
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Entrance 



Fuel region composition: 

Uranium nuclei (U + , U ++ ), percent 36 

Fiydrogen nuclei (H, H + ), percent 24 

Electrons, percent 40 

Entrance: 

Rate of fuel injection, rfip, g U/sec 1.1 

Rate of propellant flow (at 300 K), ii H , g H/sec 34 

Rate of seed flow, n|j, g ““U/sec 11 

Core: 

Uranium mass, g 375 

Fiydrogen mass, g 1 

Temperature, K 18 900 

Average dwell time, 0, sec 330 

Cavity: 

Hydrogen mass, g 172 

Uranium-238 mass, g 59 

Average temperature, K 2500 

Average dwell time, 0, sec 55 

Average exit temperature, K 3900 


(c) Average conditions in test cavity. 
Figure 1. - Concluded. 


Brazed 
joints 

Standard fuel plate J 

Outside fuel plate J 


—0.478 

— 0.201 



8.049 


\ '-0.287 (all gaps) 

'-0. 163 

13. 98 rad. CD-11408-22 
Comb is at ends only. 


Figure 2. - MTR fuel element. Fuel volume fractions: 
aluminum, 0.415; coolant, 0.583; fuel, variable. 
(All linear dimensions are in cm. ) 
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Region 

Component 

Thickness, 

cm 

Radius, 

cm 

1 

Plasma fuel 

19.2 

19.2 

2 

Propellant hydrogen 

11.28 

30.48 

3 

Aluminum cavity liner 

• 

635 

31.115 

4 

Hydrogen plenum 

2.54 

33. 655 

5 

Inner moderator containment shell 

• 

635 

34.29 

6 

Moderator 

Variable 

Variable 

7 

Outer moderator containment shell 





8 

Driver fuel 





9 

Fuel channel containment shell 





10 

Reflector 





11 

Pressure vessel 

10. 16 




Figure 3. - Spherical calculation! model of cavity test reactor. 



z 



Figure 4. - Cylindrical calculational model of cavity test reactor. (All dimensions are in cm. I 



Figure 5. - R-8 Calculational model of cavity test reactor. 
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Reactor driver core multiplication factor, 



Figure 6 . - Effect of driver fuel loading on multiplication 
constants in cavity and driver fuel regions of cavity 
test reactor at test conditions. Calculational model is 
spherical geometry and S 4 - Pj^-19 group. 



Figure 7. - Calculated specific reactivity worths of test 
cavity fuel and reactor driver fuel in cavity test re- 
actor at test conditions. 



Figure 8 . - Effect of moderator material and thickness 
on the ratio of driver power to cavity power in a 
cavity test reactor at test conditions. Calculation 
model is spherical geometry and S 4 Pj 19 group. 



Figure 9. - Required driver power generation as function of re- 
flector material and thickness in cavity test reactor at test 
conditions. Calculational model has reference model config- 
uration (except for reflector thickness) and conditions and 
has been normalized to S 4 P 1 I 9 group with cylindrical geom- 
etry. 
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Flux, neutrons/c nr- sec 


r Aluminum shell 



(a) Radial midplane. 


r Feed hydrogen 



(b) Axial midplane. 

Figure 10. - Fast and thermal flux distributions in cavity test reactor at test conditions. 
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Flux per unit lethargy 


0. 12 eV 



Figure 11. - Flux spectra in cavity and driver fuel regions of cavity test reactor at test conditions. Calculational 
model is spherical geometry and S 4 P 1 W group. 


l.OOxlO 2 



Fractional radius of test cavity fuel region 


Figure 12. - Power density distribution in the test cavity of 
cavity test reactor at test conditions. Cavity power, 2. 7 
megawatts; fuel to cavity diameter ratio, 0.63; cavity fuel 
volume, 2. 96x10^ cubic centimeters. 
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Fractional radius, R Fractional width, Y Fractional half height (measured from midplane) 

(a) Radial traverses. (b) Azimuthal traverses. (c) Axial traverse. 
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